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Abstract 



A new and intuitive perturbative approach to time-dependent quantum mechanics problems is pre- 
sented, which is useful in situations where the evolution of the Hamiltonian is slow. The state of a system 
which starts in an instantaneous eigenstate of the initial Hamiltonian is written as a power series which 
has a straightforward diagrammatic representation. Each term of the series corresponds to a sequence of 
"adiabatic" evolutions, during which the system remains in an instantaneous eigenstate of the Hamilto- 
nian, punctuated by transitions from one state to another. The first term of this series is the standard 
adiabatic evolution, the next is the well-known first correction to it, and subsequent terms can be written 
. down essentially by inspection. Although the final result is perhaps not terribly surprising, it seems to 

be not widely known, and the interpretation is new, as far as we know. Application of the method to the 
adiabatic approximation is given, and some discussion of the validity of this approximation is presented. 



The adiabatic approximation in quantum mechanics was developed in the very early days of quantum 
mechanics 1 , and is presented in many of the classic textbooks on the subject (for one of the most detailed 
presentations, see Messiah [2]). The essential idea is very simple. A time-dependent Hamiltonian has 
instantaneous eigenstates and eigenenergies, which are the solutions to the time-independent Schroedinger 
equation using H(t), where t is viewed as a parameter rather than the time. We can visualize the evolution of 
the system in terms of a graph of the instantaneous energies and eigenstates as a function of time. Suppose 
the system is initially in an instantaneous eigenstate of the Hamiltonian. Then if the Hamiltonian evolves 
slowly enough, the system will remain in the state which evolves from the initial state. However, exactly 
what constitutes "slowly enough" remains the subject of some debate. 

In the 1980s, the adiabatic approximation came to the forefront again with the work of Berry 3 , who 
showed that the wave function of a slowly-evolving system acquires a phase factor previously thought to be 
of no physical significance. Berry's phase has been observed in a wide variety of physical contexts, and has 
been proposed in many more; see the book by Shapere and Wilczek 0] for a detailed discussion and some of 
the early experimental and theoretical investigations. 

More recently, the advent of the field of quantum computing has revived interest in the adiabatic ap- 
proximation |51lrjl[7ll51l^l ll0limil2| . One reason for this is a paradigm of quantum computation based on 
the evolution of the ground state of a Hamiltonian designed to interpolate adiabatically between a simple 
Hamiltonian with known ground state and another Hamiltonian from whose ground state the answer to an 
interesting problem can be found |131 . 

A second reason is simply that in many approaches to quantum computing qubits are manipulated by 
changing an external influence (such as a magnetic field on a spin) , hopefully resulting in adiabatic evolution 
of the qubits. 

In both these situations, it is obviously important to understand the conditions under which the evolution 
of the system can be considered adiabatic. 

In this paper, we present a new perturbative approach to time-dependent quantum mechanics problems, 
which is useful in the adiabatic limit. The state of the system is written as a power series in the number 
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of jumps undertaken by the system from one instantaneous eigenstate of the Hamiltonian to another; these 
jumps are connected by an adiabatic evolution in the instantaneous state in which the system finds itself 
between jumps. To do this, we divide the time evolution into infinitesimal increments. Within one of these 
time slices the Hamiltonian can be treated perturbatively and a perturbative expansion for the evolution 
within each time slice can be calculated. The resulting series can be visualized as a diagrammatic expansion 
which is extremely intuitive and with which higher-order terms can be written down by inspection. 

The first term of the series gives the usual adiabatic expression for the final state; the second is the 
standard correction to the adiabatic evolution (see, for example, Messiah [2]). 

A somewhat simpler derivation of the final result can in fact be given by an iterative solution of the integral 
version of the time-dependent Schroedinger equation (writing the state as a sum over the instantaneous 
eigenstates of the Hamiltonian with coefficients to be determined) , but in our opinion the derivation presented 
below is sufficiently elegant and intuitive that it is worth the slight additional complexity. 

Our result can be compared to other investigations of corrections to the adiabatic approximation. Berry 
[T5] has calculated the phase change acquired by the state in "super-adiabatic" time evolution the physical 
interpretation, thus determining corrections to the Berry phase, as an iterative sequence of approximations 
to the phase. Moody, Shapere and Wilczek have gone beyond the adiabatic approximation, calculating 
perturbatively the probability of transition to other states, resulting in a sort of depletion of the adiabatic 
state. Below, we will see that our result can be shown to agree with these. Our method, however, is quite 
different, and the physical interpretation, as far as we know, is new. 

Consider a time-dependent Hamiltonian H(t). We define the instantaneous eigenstates and eigenenergies 
\n(t)) and E n {t) as the solutions to the time- independent Schroedinger equation (parameterized by time): 

H(t) \n(t)) = E n (t) |n(t)> . (1) 

For simplicity, we assume that there are no degeneracies throughout the time evolution. The state \n(t)} is 
said to have evolved from \n(Q)). 

According to the adiabatic approximation, if the Hamiltonian evolves "slowly enough" , a system prepared 
initially in an eigenstate of H(0) will, at time T, be in the state which evolved from the initial state. Let us 
see how this happens. 

If the initial state is |-0(O)) = 1 (0) ) (not necessarily the instantaneous ground state of the system), then 

|^(T))=C7(T)|0(0)), (2) 

where the time evolution operator is 

U(T) =f exp-i [ dtH(t), (3) 



in units where fi = 1 and where T denotes time ordering. 

We divide the total time interval T into N intervals of duration e = T/N; ultimately we will take the 
limit N — > oo, e — * with T fixed. Define tj = je, the beginning of the jth time interval. Then we can write 

U (T) = Un-iUn-2 ■■■UiUo, (4) 

where Uj is the evolution operator taking the system from tj to tj+\\ 

Uj=f exp-i ( 3+ dtH(t). (5) 

Let us find an expression for the state after one time interval, \tj}{t{)) = Uo \4>(0)) — Uq |0(0)). It is useful 
to expand this state in eigenstates of H(t±): 

|V(ti)) = ^TM^MMIt/oloto)) 

n i 

/■ti 

= V \m(ti)) (m(ti)|Texp-« / dt H(t) |0(0)) . (6) 
™ Jo 
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So far, no approximation has been made. We will now do so. Writing H(t) = J2 n l n W) E n (t) (n(t)\, we 
expand each of these factors about t = 0: 

|n(t)> = |n(0)) + t|n(0)> + ... 

E n (t) = E n {0) + tE n {0) + ■ ■ ■ (7) 
(n(t)| = (n(0)| + t(n(0)| + -». 

Since we are interested in times less than e, which eventually will tend to zero, as long as all derivatives 
are finite the series will converge for sufficiently small e. (Note that we have not as yet assumed adiabatic 
evolution.) Then 

H 



(*) = £ K"))^(0)('»(0)l+<(WO))^(0)(n(0)l 

n ^ 

+ |n(0)> £7„(0) (n(0)| + |n(0)> E n (0) (n(0)|) + 0(i 2 ) j (8) 



In fact, as we shall see shortly, only the first term is necessary. 
We can now write down an approximation of Uq in (JBJ: 



U = l-i I dtH(t)+0(e z ) 
Jo 

= l-teJ2n0))E n {0){n(0)\+O(e 2 ). (9) 



Note that the 0(t) terms in (jSJ give rise to 0(e 2 ) terms in 10. 
We also expand {m(ti)\: 

(m(tx)\ = (m(0)| + £ (to(0)| + 0(e 2 ). (10) 
The matrix element in © can now be written 

(m{h)\U o \0(0)) - ((m(0)| + e(m(0)|) ^1 - iej^ \n(fl)) En(Q) (n(0|^ |0(0)> + 0(e 2 ) 

= <5 m0 - ieE n (0)d m0 + e (m(0)| 0(0)} + (9(e 2 ) 

'( )<5 m o + e(m(0)|0(0))+0(e 2 ). (11) 



e 



-ieE, 



Substituting this in ©, 

m,)) = e~^(°) |0(t x )) + eJ2 (m(0)| 0(0)} |m(tx)) + 0(e 2 ) (12) 

m 

Now, we are free to choose the phase of the instantaneous eigenstates as we please; we can use this freedom 
to set (n(t) \ n(t)) = for all n, 1 so that with this choice the sum over m in l|12ll does not include m = 0. 

We have succeeded in writing \tp(ti)} in terms of instantaneous eigenstates of H(t\) with the coefficients 
given by power series in the infinitessimal parameter e. The dominant term (as long as the derivatives of H 
are finite), not surprisingly, is |0(ti)) - not because of any assumed adiabaticity, but simply because we have 
not given the system enough time to have much of a chance of changing states. 

We can now apply this procedure to the second time interval, and so on. We can look at this as sub- 
stituting (@J into J2J> inserting 1 = J2m j+1 l m j+i(^'+i)) ( m j+i(^'+i)l to tne ^ 01 eacn Uj, and expanding 

1 Locally, this is true; that we must be careful globally was discovered by Berry the integral of this factor over a closed 
loop in the parameter space of the Hamiltonian is none other than Berry's phase. 
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everything in sight, rather reminiscent of derivations of the path integral form of the propagator in quan- 
tum mechanics (where, of course, position eigenstates are inserted). This results in the following fairly 
cumbersome expression: 



N-l 

IV>CT)> = £ E ••■EEl m ^( T )> II (W^e- i ^( t ^+e<m i+1 fe)|m i (t j )) + 0(e 2 )), (13) 

rriN ttim-1 ma m± j—0 

where too is not summed over: too = 0. For small e, in each of the factors in parentheses the first term 
is obviously dominant; nonetheless, in the product of these factors there are so many more contributions 
coming from the e term that it is just as important when e — > 0. We will show, however, that the e 2 term is 
indeed negligible. 

Let us now dissect by keeping track of the number of e terms we include in the product, the remaining 
factors being the first term. To begin, we take the first term from every factor in the product. We obtain 
what will be referred to as {^(T)}^ : 

N-i 

iw) (0) = E E ••■EEi™^)> II w.^-**"^- ( 14 ) 

rnjs; mjv-1 m2 mi j—0 

The (S-functions eliminate all the sums, and in the limit TV — > oo we get 

|-0(T)) (O) = e' i Io dtE ^ t) \0(T)). (15) 

This term is the usual adiabatic result: the state is the instantaneous eigenstate which evolved from the 
initial state, multiplied by the expected dynamical phase factor. 

The next contribution to \ip(T)), referred to as \ip(T))^ , is obtained by taking one e (to| to) term and 
otherwise taking the first term from the product in i|13|) . We must sum over which factor contributes the 
e (to| to) term (the kth factor, say): 

AT-l / 

l^( T )) (1) = E E E •••EEl m ^ T )> [T[Sm j+1 , mj e- UE ^ } e(m k+1 (t k )\m k (t k )). (16) 

k=0 mm mjv-i m 2 mi 

Only the sum over m k +\ survives; earlier TOfc's are set to zero while later ones are set to m k +i- Renaming 
tnk+i ~^ m and t k — > t\ and taking the limit N — * oo, we obtain 

MT)) ll) = E M T )> l T dt ie -^ dtE - {t) (m(t 1 )|0(t 1 )>e- < r**«. (17) 

This contribution represents an adiabatic evolution in the initial state (more precisely, the state evolving 
adiabatically from the initial state) until time t\, a transition — » to at that time, followed by another 
adiabatic evolution in the state \m). So \ip{T)y' is the "one-jump" contribution to \ip(T)) (whence the 
superscript). 

The next contribution to \ip(T)), taking two e(rh\m) terms, is the "two-jump" contribution; the result 
can be written by inspection: 

|^(T)} (2) = EEM T )> f 2 dt ie - l ^ dtEm{t) (TO(t 2 )|n(t 2 )) 



m^n n^O 

e 



l ^ dtE " {t) (n(t 1 )\0(t 1 )}e- t C dtEoit \ (18) 



The physical interpretation is clear: there are transitions at times t\ and i 2 > ty, otherwise, the evolution 
is adiabatic. The generalization to higher terms is obvious. 
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It remains to demonstrate that the 0(e 2 ) terms in Ill-iH are unimportant. This can be shown by straight- 
forward power counting. For instance, if we consider taking one 0(e 2 ) term and otherwise taking the first 
term, we get an expression of the form 

e2 EEM T )> cx P \- ie E E mitj)) X k exp (-ief^Eoitj)) , (19) 

k rn \ j=k+l J \ j=0 ) 

where Xf. contains various matrix elements, energies, etc. We can convert e^2 k and e t° integrals, leaving 
one factor e left over, indicating that this term is vanishingly small as e — > 0. 
Putting all the terms together, we get an expansion for \ip(T)): 

|V>(T)) = |^(T)) (0) + |^(T)) (1) + |VXT)) (2) + • • • . (20) 

We can give a diagrammatic representation for each term in the series. With time going to the right, 

n 



_L 







+ ••• (21) 



where each horizontal line represents a factor exp— i JdtE n (t), each jump n — > m represents a factor 
(m(t) \ n(t)), and we add a ket corresponding to the final state at time T (|0(T)) in the first diagram, \m(T)) 
in the rest). Finally, the times of the jumps must be integrated over (respecting the ordering) and the states 
m, n, etc. must be summed over. 

The expansion ill'OII can also be written in the following compact (albeit not terribly useful) form: 

|V(T)) - M T )> u mo(X), (22) 

m 

where U is the time evolution operator written in the "moving" basis given by the instantaneous eigenstates. 
This is 

U mn {T) = e- 1 $o Em(t)dt (f expx) , (23) 
where ^ 

dt (m{t)\n{t))e- l ^ (E - {t) - Em(t))dt . (24) 



X 



inn 







In this form, our result can be seen to agree with that of |16|. and with |15j for the no-transition case. 

We can easily calculate the amplitude to end up in the state |m(T)), (m(T)\ ip(T)); for this we can use 
the same diagrammatic expansion but we do not sum over the final state and do not include the final ket in 
the analytic expression. For example, the amplitude to end up in the state |0(T)} (that is, to have made no 
net jump) is 



<0(T)|V(T)) = 







-if o T d tE o{ t)( 1 + y2 f T dt2 f 2 dtl {b{t 2 )\n{t 2 )) 



-i J' 2 dt(E n (t)-E (t)) 

xe J *i (n(ti)|0(*i))H } (25) 
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while the amplitude to have jumped to the state \m(T)} is 



<m(T)|V(T)> = 



m 

o r ~ 



n 

1 m 



+ 



ti t\ t2 

«j^**(t}| f T dti -is:^iE m{ t)-E 0{t)) {fh(ti) | 0(tl)) 







+ 2^ / "*2 / dt\& Jt 2 (m(t 2 )| n(i 2 )) 

/n v/0 JO 

ri7=0,m 

-t f* 2 dt(B„(t)-B (t)) 1 

xe Jf i (n(ti)|0(t 1 )}H ^ (26) 

The expansion H2Q(I for |-0(T)) is obviously an expansion in powers of the number of jumps, or transitions, 
the system undergoes. Under what circumstances is this expansion useful; that is, under what circumstances 
will it converge? Intuitively, we expect this to occur when the Hamiltonian is slowly-varying since then 
transitions are expected to be rare. In Eqs. I|15I17I18[) the presence of increasing powers of matrix elements 
(rh\ n) suggests that indeed as the Hamiltonian varies more and more slowly, the instantaneous eigenstates 
will change more and more slowly, and the expansion should indeed get more and more convergent. 

However, it should be noted that the smallness of the matrix elements in (|26|) may be cancelled by the 
largeness of the interval over which the time is integrated. It is in fact the relatively slow variation of the 
matrix elements on the time scale related to the Bohr frequency that is the key to improved convergence of 
(I2t)(l . and to the validity of the adiabatic approximation. 

Before examining the adiabatic case, let us look at the case of non-adiabatic evolution. A trivial obser- 
vation is that for a given Hamiltonian H(t), then for T sufficiently small, integrals such as those in (|17I18|) 
go to zero, and the system is to a good approximation in its initial state, for the simple reason that we have 
not given it enough time to make a transition. The important question is: How short must T be for this to 
be so? If the scale on which the Hamiltonian varies is t, then generically | (rh\ n)\ ~ (t) _1 , and dimensional 
analysis indicates that 1 1 \ip (t)) (k) \ | ~ {T/r) k . Thus, successive terms in (|20|l are smaller if T < r, and 
the series will converge (at least naively; a rigorous statement about convergence would of course require a 
detailed analysis). We can conclude that the initial state will be depleted on a time scale T ~ r, which is 
as one might expect: the time scale of depletion of the initial state is none other than the time scale of the 
Hamiltonian. 

So far, we have made no use of the other time scale available, namely, the inverse of the characteristic 
energy scale of the problem. This is of course the key to the adiabatic approximation. We have three time 
scales available: the inverse of the characteristic energy scale (E, say) of the problem, the rate of change of 
the Hamiltonian r, and the duration of the process T. It is common to consider the latter two time scales 
to be equal, but for generality we will not do so here. 

To more easily keep track of the relative size of the various terms, it is useful to make a change of notation 
to dimensionless quantities. The Hamiltonian is characterized by a time scale r and an energy scale E, so 
we define a dimensionless time s — t/r, and a dimensionless Hamiltonian and its energies and eigenstates 
which depend smoothly upon s; these are h(s) = H(i)/E, e n (s) = E n (s)/E and \n{s)) = \n(t)). We make 
the following assumptions: 

1. The state |0(s)) has zero energy, eo(s) = 0. (This assumption is for simplicity only and can easily be 
removed with only minor modifications to the analysis.) 

2. The minimum energy gap to any other state is unity: |e n (s)| > 1 V n =/= 0. 

3. Otherwise, h(s) is "generic": derivatives of the energies e n (s) are of order 1 (for n 0, of course), 
and derivatives of the states \n(s)) have norm of order 1. (This is simply the statement that H(t) is 
characterized by time scale r, recast in terms of dimensionless quantities.) 
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The duration of the process is T, so we define a dimensionless duration S — T/t. We define a parameter 
A by A = Et\ when A ^> 1, the Hamiltonian can be said to evolve slowly, and the evolution of the system 
is expected to be adiabatic. We also define a parameter 7 by 7 = max s;m ^o I (™'( s )l 0(s))|, where the prime 
denotes differentiation with respect to the argument; 7 is of order unity. 

We can now rewrite the expansion l|2U|) in terms of dimensionless quantities. For instance, the first 
correction is, from lfT7|) . 

|^(S)>« - V \m(S)) f S d Sl e~ lX ^ dsem(s) (m'( Sl )|0( Sl )>. (27) 

Since the zeroth term \tp(S))^ is of norm unity, a first indication of the convergence of the series, and the 
validity of the adiabatic approximation, comes from the norm squared of the first-order correction: 

W{fcS)\#(S))W = £|A«(S)| 2 , (28) 

where ^ 

A$(S) = (m(S)\i,(S))M = f d Sl e~ lX ^ dsemis) (m'( Sl )\6( Sl )). (29) 

Jo 

Note that the smallness of \A m ^\ is not due to the smallness of the matrix element (which is of order one). 
To see that \Am\ goes to zero in the limit A — > 00, we note that since e m and its derivatives are of order 1, 
the exponential oscillates rapidly, on a (dimensionless) time scale (A) -1 . The matrix clement, in contrast, 
varies on time scale 1, so that within one period of the exponential the matrix element is approximately 
constant, and the exponential integrates to zero. Corrections to this give a result of order 1/A. 
As a special case, if e m (s) and the matrix element are both constants (equal to one, say), then 

1 1 — e~ iXS I 

= J ~ x L > (30) 

which clearly goes to zero as 1/A. 

To analyze Am in the general case, we define f m (s) = /„ s ds'e rn (s') and g m {s) = (m'(s)\ 6(s)). Then 

A (i) = e -ix fm (s) [ S dse ^U(s) gm{s y ( 31) 



The assumptions stated above tell us that f m , g m and their derivatives are of order 1, and that \f m \ > 1. It 
is then easy to show that if S ~ 1 and A>1, 

A® = ^ (HS -e-^ 9 -^§-) +0(A- 2 ). (32) 

This clearly demonstrates that the dominant contribution to | is the first term when A > 1; rewriting 
it in terms of dimensionful quantities, 

^^^^-e-'SI^^mmL^^). (33) 
iE m (T) iE m (0) 

Interestingly, the dominant correction depends only on the initial and final values of the matrix element; 
in particular, in the special case where initially and finally the instantaneous eigenstates are constant, the 
dominant correction is zero and the adiabatic approximation should be more robust, being valid up to 
corrections of order l/(ET) 2 rather than l/(ET). 



7 



Assuming this special case does not apply to the system under investigation, we can now give a bound for 
\A^m\ in terms of quantities defined earlier. The matrix elements satisfy \g m \ > 7, while the dimensionless 
energies e m are at least unity. So 



Thus, we would conclude on the basis of this first-order analysis that the adiabatic approximation is valid 
if the right hand side of 1)34(1 is small; this is essentially the "standard" condition of validity of the adiabatic 
approximation which has been re-examined recently (51 1?>1 171 l8l l9l 1101 1111 112) . This leads us to the obvious 
question: Docs the first-order analysis extend to higher orders; that is to say, is \Am | ~ {Er)~ k 7 As 
we will see, the answer is no. In a way, this is not so surprising: if indeed Am behaved in this way, we 
would conclude that a slowly-evolving Hamiltonian will never give rise to a transition with any appreciable 
probability. Surely this is not so; for fixed rate of change of the Hamiltonian, the system ought to eventually 
escape from the instantaneous eigenstate related to the initial state. 

The fundamental reason why the first-order analysis does not extend directly to higher orders is most 
easily seen by looking at A^ , the second term in i|25|) . We can write the dimensionless version of this as 



From (|32[1 . we see that An contains a non-oscillating piece, and the argument used to show the smallncss 
of l(2"§|) does not apply to (J35J; indeed, 



Rather than being a factor (Et)^ 1 times the first-order term, it is a factor of order T/t times the first-order 
term. Taken at face value, we would conclude that the expansion no longer converges, and presumably the 
initial state is depleted, on a time scale T ~ r: we would be no better off than the non-adiabatic evolution. 
However, the situation is not quite as dire as this: higher-order terms can be shown to get smaller up until a 
time of order T ~ Et 2 . Thus, the depletion time scale appears to be longer than the non-adiabatic case by 
a factor Et 3> 1 (the inequality being true for a slowly- varying Hamiltonian) , so indeed some time is gained 
in the adiabatic limit. This will be analyzed in more detail in a future publication |17| . 

As a final remark, we suspect that one factor giving rise to the apparent inadequacy of the "standard" 
condition of validity of the adiabatic approximation lies in the fact that the examples given in some papers 
questioning this condition may not fit into the type of problem discussed in this paper (and in at least some 
of the standard discussions of the adiabatic method, such as that of Messiah [2]), that is to say, they are not 
described by a Hamiltonian which can be written H(t) = Eh(t/r), where h(s) is "generic" as defined above. 
More precisely, a Hamiltonian containing a part which is rapidly oscillating but small in amplitude can 
satisfy the "standard" condition, yet the evolution might not necessarily be adiabatic; such a Hamiltonian 
is not "generic" as defined above. This is currently under investigation |17|. 

Note added: After this paper was submitted, a preprint by T. Vertesi and R. Englman has appeared, 
which studies possible failures of the standard condition of validity of the adiabatic approximation from a 
perturbative point of view |18| . 

We thank A. A. Methot, M.B. Paranjape, and B. Sanders for useful conversations. This work was funded 
in part by the National Science and Engineering Research Council. 
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